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ABSTRACT 


This  document  describes  a  FORTRAN  IV  subroutine  named  CAVAP 
(an  acronym  for  conCAVe  Assignment  Problem)  which  uses  a  branch 
and  bound  method  to  solve  a  variant  of  the  assignment  problem  in 
which  the  objective  to  be  minimized  is  concave.  It  is  assumed 
that  m  tasks  are  to  be  performed  and  n  different  kinds  of 
machines  have  been  defined  for  possible  use  in  performing  them. 

In  general  a  task  may  be  performed  by  several  machines  of  one 
type,  or  a  different  number  of  machines  of  another  type,  or 
possibly  only  by  certain  combinations  of  machines.  It  is  assumed 
that  total  cost,  which  is  to  be  minimized,  is  the  sum  of  costs 
for  each  kind  of  machine,  each  of  which  is  in  turn  a  concave 

function  of  the  number  of  machines  of  that  type  required.  .  A 

y 

detailed  mathematical  treatment  of  the  problem  and  the  solution 
method  has  been  given  elsewhere _  •  CAVAP  has  been  used  to  solve  a 
space  fleet  selection  problem  with  20  space  missions  (tasks), 

40  booster  and  spacecraft  components  (machines),  and  a  number  of 
alternative  launch  vehicle  configurations  and  combinations  for 

each  mission.  Running  time  was  under  20  minutes  on  an  IBM  360 

/\ 

model  44  with  one  disk. 


1.0  GENERAL  INFORMATION 

This  document  describes  a  Fortran  IV  subroutine  named  CAVAP 
which  solves  a  variant  of  the  of  the  so-called  assignment  problem 
in  which  the  objective  to  be  minimized  is  concave.  CAVAP  is  an 
acronym  for  conCAVe  Assignment  Problem.  The  method  of  solution  is 
by  brancn  and  bounds.  A  mathematical  treatment  of  the  problem  and 
solution  method  is  given  in  [2].  For  an  example  of  the  use  of  this 
subroutine,  see  the  description  in  [1]  of  the  routine  for  computing 
a  minimum  cost  mixed  space  fleet  from  given  vehicle  designs  when 
the  fixed  costs  and  learning  curves  give  a  concave  cost  vs.  quantity 
relationship.  CAVAP  has  been  run  as  a  subroutine  to  the  program  in 
[1]  on  an  IBM  360  model  44  computer  equipped  with  one  disk.  All 
limitation  and  running  time  data  presented  are  taken  from  this 
context. 

1.1  Purpose 

1.1.1  Mathematical  statement  of  the  problem 

Suppose  m  tasks  are  to  be  performed  and  n  different  kinds 
of  machines  have  been  defined  for  possible  use  in  performing  them.  In 
general  a  task  may  be  performed  by  several  machines  of  one  type,  or  a 
different  number  of  machines  of  another  type,  or  possibly  only  by 
certain  combinations  of  machines.  For  each  i,  1  <_  i  m,  it  must 
be  possible  to  write  out  all  the  alternatives  to  be  considered  for 
performing  task  number  i.  Let  k(i)  denote  the  number  of  such 
alternatives.  For  each  value  of  h,  1  <  h  <_k(i),  the  h^  alternative 
for  task  i  can  be  given  as  a  vector 
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where  each  is  the  number  (positive  or  zero)  of  machines  of 

type  j  required  in  alternative  h  for  performing  task  i.  There 
is  no  requirement  that  the  alternatives  for  different  tasks  are  in 
any  way  related.  For  the  application  given  in  [1],  the  tasks  are 
space  missions,  perhaps  20  in  number;  the  machine  types  are  booster 
and  spacecraft  components,  perhaps  40  in  number;  and  each  alternative 
consists  of  the  numbers  of  components  of  different  types  needed  to 
build  the  number  of  launch  vehicles  of  a  particular  configuration 
required  to  perform  one  mission. 

A  feasible  (not  necessarily  optimal)  solution  of  a  problem 
consists  in  selecting  exactly  one  alternative  for  each  task,  i.e. 
numbers  h(i).  It  is  assumed  that  a  machine  to  be  used  on  one  task 
cannot  also  be  used  on  another,  but  different  machines  of  the  same 
type  may  be.  Thus  the  number  x(j)  of  machines  of  type  j  required 
by  a  particular  feasible  solution  is  the  sum  of  the  numbers  of  machines 
of  type  j  included  in  the  selected  alternatives,  i.e. 


x 


xih(i) 

xj 


Next  suppose  the  total  cost  of  design,  production,  etc.,  of  x(j) 
machines  of  type  j  is  given  by  a  cost  function  4>(j,x(j)).  This 
implies  that  the  cost  of  x(j)  machines  of  type  j  is  completely 
independent  of  the  numbers  or  kinds  of  other  machines  produced.  Thus 
in  the  problem  given  in  [1],  a  complete  vehicle  cannot  be  considered 


cost 


a  machine  if  for  example  it  has  a  stage  in  common  with  another  vehicle. 

It  is  also  assumed  that  for  each  j  the  graph  of  $(j,x(j))  is  concave, 
that  is,  the  cost  of  an  additional  machine  of  type  j  decreases  as  the 
total  number  of  machines  of  type  j  increases.  Figure  1  illustrates 
acceptable  and  unacceptable  functions  <(>.  In  most  practical  applications 
it  is  to  be  expected  that  $  will  be  zero  for  x  *  0  and  positive 
increasing  for  x  >  0  as  in  Figure  lb,  but  this  is  not  essential  to  the 
method  of  solution.  Since  no  alternative  can  specify  a  negative  number 
of  machines  of  any  type,  it  is  immaterial  what  value,  if  any,  <C  gives 
when  x(j)  is  negative. 


Figure  1 

The  total  cost  of  performing  all  tasks  using  a  particular  feasible 
solution  is 


Z 


j-i 


The  problem  solved  by  CAVAP  is  that  of  selecting  the  feasible  solution 
(i.e.  an  alternative  for  each  task  or  h(i)  for  each  i)  which  will 


minimize  Z. 


1.1.2  Need  for  a  special  algorithm 

The  method  of  complete  enumeration  is  an  obvious  and  very  simple 
one  for  determining  the  least  cost  solution  of  any  assignment  problem. 
This  is  the  method  which  computes  the  cost  of  each  of  the  possible 
combinations  of  alternatives  and  retains  the  combination  with  the 
least  cost.  This  method  requires  the  examination  of 

s  *  k(l) *k (2)  • .  . .  *k(m) 

different  feasible  solutions.  For  one  of  the  problem?,  of  the  type 
described  in  [1],  m  ■  20,  n  *  38,  and 

s  -  5*l*5*3*12*5*ll*2*6*7*5*6*15‘6-15*6*6*3-7*7  -  101:>. 

Even  if  the  cost  of  each  feasible  solution  could  be  computed  in  1 
microsecond,  it  would  take  more  than  30  years  to  examine  all  lO1^ 
of  them.  By  comparison  the  same  problem  was  solved  using  CAVAP  in 
approximately  20  minutes  on  an  IBM  360  model  44. 

1.2  Limitations 

The  method  of  solution  used  by  CAVAP  belongs  to  a  broad  class 
called  "branch  and  bound"  methods.  As  a  class  these  methods  suffer 
two  inherent  limitations.  First,  as  the  size  of  the  problem  increases, 
the  amount  of  computation  time  required  increases  exponentially.  This 
means  that  there  is  a  fairly  well  defined  bound  on  problem  size  above 
which  computation  time  is  impractically  large.  The  better  the  procedure, 
the  higher  this  bound.  Second,  branch  and  bound  methods  make  use  of  a 
list  of  partial  solutions.  In  order  to  save  recomputing  time,  the 


amount  of  information  saved  with  each  item  of  the  list  is  large. 

Thus  auxiliary  memory  is  usually  used.  Since  the  length  of  the  list 
grows  exponentially  with  the  size  of  the  problem,  there  is  also  a 
rough  limit  on  problem  size  above  which  the  list  is  too  long  even  for 
auxiliary  memory. 

For  the  application  discussed  in  [1]  it  was  found  that  the 
limitation  on  problem  size  was  the  size  of  the  scratch  area  on  the 
only  available  disk,  namely,  500  blocks  of  180  words  each.  The 
largest  problems  which  did  not  exceed  this  storage  requirement  ran 
about  20  minutes.  For  several  scratch  disks  the  limitation  would 
become  computation  time. 

Presumably  with  large  problems  round-off  errors  could  become 
serious.  In  many  computer  programs  the  effect  of  round-off  errors 
is  to  cause  a  fairly  steady  deterioration  of  accuracy  as  the  amount 
of  computation  increases.  The  nature  of  the  algorithm  used  by  CAVAP 
is  such  that  with  increasing  size  of  problem  there  will  be  little 
effect  on  the  precision  of  answers  until  at  some  point  round-off 
errors  cause  the  process  to  go  wild.  It  is  unlikely  that  this  would 
result  in  undetected  errors.  The  most  likely  result  would  be  an 
immediate  error  exit  or  an  endless  loop  which  would  eventually  fill 
up  available  disk  storage  and  result  in  an  error  exit.  Results  with 
the  problems  of  [1]  suggest  that  the  limit  due  to  running  time  will 
be  met  before  round-off  beciines  a  problem. 

As  written,  CAVAP  will  handle  in  theory  up  to  44  machines,  up 
to  100  tasks,  and  up  to  500  auxiliary  memory  items.  These  limits  are 
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easily  alterable.  The  only  real  limitations  are  the  amount  of  running 
time  and  auxiliary  memory  available. 

2.0  PROCEDURE 

The  naming  of  variables  in  the  subroutine  follows  closely  the 
notations  found  in  [2]  and  in  paragraph  (1.1.1). 

The  method  of  solution  is  given  in  considerable  detail  in  [2]. 

The  deviations  are  minor  and  of  a  programming  nature.  The  subroutine 
contains  an  initial  section  in  which  the  input  parameters  and  data  are 
subjected  to  a  complete  check  to  Insure  that  they  meet  the  requirements 
set  fourth  in  paragraph  (3.0)  below. 

In  addition  to  the  normal  return  there  are  two  abnormal  returns 
from  the  subroutine.  If  the  normal  return  occurs,  a  feasible  solution 
has  been  found  which  has  a  cost  which  differs  from  the  minimum  cost  in 
about  the  5th  decimal  place.  Abnormal  return  1  occurs  if  the  input 
data  is  Improperly  given.  If  the  program  must  be  abandoned  for  any 
reason  before  finding  an  optimum  solution,  a  feasible  solution  and  its 
cost  are  given.  It  may  very  well  happen  that  the  feasible  solution  in 
this  case  is  actually  optimum,  but  the  subroutine  has  not  yet  succeeded 
in  proving  this  fact.  Three  likely  causes  of  the  subroutine  being 
abandoned  are:  exhausting  the  disk  storage,  finding  that  one  of  the 
cost  functions  given  by  PHI  is  not  concave,  or  the  accumulation  of 


round-off  errors. 


-7- 


3.0  INPUT  AND  OUTPUT  DESCRIPTION 
The  calling  sequence  is 


CALL  CAVAP(M,N,PHI,JIN,XIN,ZEEST,HIBEST,XBEST,NITEf  ,MAXITM,&a.6b) 

» 


where 


Name 


M 

N 

PHI 


JIN 

XIN 


ZBEST 


HIBEST 


XBEST 


Type 
Integer 
or  Real 


I 

I 

R 


I 

R 

R 

I 


R 


Description 


[  Input) 

number  of  tasks.  MLSt  be  between  1  and  100. 
maber  of  machines.  Must  be  between  1  and  44. 
transmitted  name  of  a  function  subroutine 
PHI(J,X)  which  returns  the  cost  ot  X 
machines  of  type  J.  This  function  »■  st 
be  concave  over  the  range  of  values  of  X 
for  the  given  problem.  The  function  need 
not  be  defined  for  negative  X. 

input  arrays  described  below. 

(Output) 

cost  of  the  best  feasible  solution  found.  Set 
to  Infinity  if  the  input  data  is  improper, 
an  array  whose  1th  element  is  the  number  h 
of  the  alternative  selected  for  task  1  in 
the  best  feasible  solution  found.  This  array 
should  be  of  length  at  least  100  M. 
an  array  whot»e  Jth  clement  is  the  number  of 


machines  of  type  j  required  in  the  best 
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feasible  solution  found.  This  array 
must  be  of  length  at  least  N. 

NITER  I  number  of  iterations  performed. 

MAXITM  I  maximum  number  of  items  on  disk  at  any  time 

during  solution  of  current  problem. 

(Error  returns) 

a  statement  number  of  return  if  input  data  is 

improper. 

b  statement  number  of  return  if  solution  is 

abandoned. 

The  input  arrays  JIN  and  XIN  contain  the  descriptions  of  the  various 
alternatives  for  each  of  the  tasks  arranged  in  a  compressed  form.  The 
first  element  of  JIN  must  be  -1.  The  second  element  of  JIN  must  be  0. 
There  must  follow  a  string  of  from  1  to  N  distinct  positive  integers 
less  than  or  equal  to  N.  These  are  the  numbers  of  the  machines  used  in 
the  first  alternative  for  the  first  task.  In  the  corresponding  locations 
in  XIN  must  be  the  "numbers"  of  these  machines  required.  These  numbers 
must  be  floating  positive  numbers;  they  need  not  be  floating  integers. 

The  first  location  of  JIN  following  this  string  must  contain  0  or  -1. 

If  the  entry  is  0,  another  string  describing  a  second  alternative  must 
follow.  Any  number  of  alternatives  may  follow,  each  one  preceded  by  a 
0  in  JIN.  If  the  entry  following  the  description  of  an  alternative  is 
-1  the  description  of  the  alternatives  for  the  next  task  must  follow. 
Again  the  first  alternative  must  be  preceded  by  a  0  entry  in  JIN.  The 
first  entry  in  JIN  following  the  last  alternative  of  the  last  task  must 
be  -2.  Each  task  must  have  at  least  one  alternative  and  the  number  of 
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tasks  described  (number  of  -1  entries  in  JIN)  must  equal  M.  An 
example  of  an  acceptable  set  of  values  for  JIN  and  XIN  involving  2  tasks 
with  1  and  2  alternatives  respectively  is  shown  below.  The  value  of  N 
is  A.  The  values  marked  *  are  immaterial. 


JIN  XIN 

-1  * 

0  * 

1  2.0 
A  1.5 
-1  * 

0  * 

2  3.0 

0  * 

3  1.2 

2  .05 

A  5.0 

-2  * 


oooooooooo 
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FORTRAN  IV  MODEL  44  PS  VERSION  3.  LEVEL  2 

SUBROUTINE  CAVAP  I M, N, PHI , J IN, X  IN, ZBEST , H IBEST, XBEST, 

1  NITER,  MAX  I TM, *, * ) 

FOR  A  DESCRIPTION  OF  THIS  SUBROUTINE.  SEE 

•CAVAP.  A  FORTRAN  SUBROUTINE  FOR  SOLVING  AN  ASSIGNMENT 
PROBLEM  WITH  CONCAVE  OBJECTIVE' ,  BY  D.  WALKUP. 
MATHEMATICAL  NOTE  NO.  529,  BOEING  SCIENTIFIC  RESEARCH 
LABORATORIES.  SEPTEMBER  1967. 

FOR  A  MATHEMATICAL  TREATMENT  OF  THE  PROBLEM  ANO  THE  THEORY  OF  THE 
ALGORITHM  USED,  SEE 

•ON  A  BRANCH- AND- BOUND  METHOD  FOR  SEPARABLE  CONCAVE 
PROGRAMS*,  BY  D.  WALKUP,  MATHEMATICAL  NOTE  NO.  527, 
BOEING  SCIENTIFIC  RESEARCH  LABORATORIES,  SEPTEMBER  1967. 
IMPLICIT  INTEGER  IH-N) 

EXTERNAL  PHI 

DIMENSION  J1N(1),XIN(1),HIBESTI1) , XBEST ( 1 ) 

DIMENSION  HIL INI  100) , XL  INI 44) ,2  1 44 1 , ZL INI T ( 5001 ,LOC ITMI 500 1 
DIMENSION  INFO!  180 ) , P I  44 ) ,Q( 44 ) , C I  44 ) ,01 44 ) 

EQUIVALENCE  I  INF 0(1), PI  IP, I  INFO  t  45)  .  01  1 )  )  • 

1  I  INFO!  89)  ,CU) ),  IINFOI 1 33)  ,01 1 ) )  ,  IINFOI 177 ) ,  DSUM I , 

2  I  INFO! 1 7B ) , JBK  ) , I  INFO! 179) , XBK ) . I  INFO! 180) ,ZBK) 

DIMENSION  XIH(44)*X1MAX(44) 

EQUIVALENCE  <P,X1H),(XLIN,XIMAX) 

DATA  FINF/Z7FFFFFFF/ 

DEFINE  FILE  1 1 1300,90 , U.  IX ) 

XEPS  >  .00001 
OMEPS  -  .99999 
MMAX  >  100 
NMAX  -44 
ITEMAX  «  500 
LINFO  -  180 
NITER  «  0 
NITEMS  «  0 
MAXI TM  «  0 
ZBEST  •  FINF 
L  «  0 

IFI (M.LE.O).OR. (M. GT. MMAX). OR. ( N.LE.O ) .OR . ( N.GT. NMAX) )  GO  TO  900 
C  CHECK  INPUT  LISTS  JIN  ANO  XIN  AND  FIND  UPPER  BOUNDS  Of J)  ON 

C  THE  SOLUTION. 

DO  105  J  «  1 , N 
105  Q I J  )  >  0. 

I  •  0 
L  »  1 

IF!  JINIU.NE.-l)  GO  TO  900 
1101  *  I  ♦  1 

IF(I.GT.M)  GO  TO  900 
00  115  J  «  1 . N 
115  X  I  MAX  I J )  •  0. 

H  -  0 
L  -  L  ♦  l 


n  n  o  n  o  oooo 
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IM  JIN(L).NE.O)  GO  TO  900 
120  H  «  H  ♦  l 

DO  125  J  *  1  »N 
125  XIH1JI  =  0. 

IF(JlN(L*l l.LE.OI  GO  TO  900 
130  L  «  L  ♦  1 
JJ  *  J 1 N ( L  > 

IF(JJ.LE.O)  GO  TO  140 

I F ( (JJ.GT.N).OR.(XIN(L).LE.O. )  • OR .(XIH(JJ).GT.O.))  GO  TO  900 
XIH(JJ)  *  X  I N ( L ) 

GO  TO  130 

140  DO  141  J  «  1  ( N 

IF(XIH(J)  .LE.XIMAX  ( J )  )  GO  TO  141 
X  I  MAX ( J )  »  X  I H( J I 

141  CONTINUE 

145  IF(JJ.EQ.O)  GO  TO  120 
DO  150  J  *  1  * N 
150  OIJ)  «  Q ( J )  ♦  X 1  MAX  (  J ) 

IF< JJ. 60.-1  )  GO  TO  110 
IF(I.NE.M)  GO  TO  900 

CHECKING  OF  INPUT  LISTS  JIN  AND  XIN  COMPLETE.  USE  OF  XIH  AND 
X  I  MAX  FINISHED. 

STEP  1.  FIND  FIRST  LINEARIZED  SOLUTION. 

COMPUTE  INITIAL  LINEARIZED  OBJECTIVE  FUNCTIONS. 

DSUM  *  0. 

DO  1 60  J  *  1 1 N 
PIJI  *  0. 

ZP  *  PHKJ.O.) 

D(J)  *  ZP 

DSUM  =«  DSUM  ♦  ZP 

C(J)  «  0. 

IF (Q( J) .EQ.O. )  GO  TO  160 
IQ  *  PHllJ.QiJI) 

C(J)  «  (ZQ-ZPJ/QIJI 
160  CONTINUE 

INITIALIZE  DISK  LOCATION  LIST 
DO  170  L  *  1 1  I TEMAX 
170  L0C1TM1L)  «  2*L-l 

SOLVE  INITIAL  L (NEAR  I Z  EO  PROGRAM  AS  A  SPECIAL  CASE  OF  STEP  4. 
NCASE  *  2 
GO  TO  240 

STEP  2.  SELECT  ITEM  FROM  DISK. 

(REMOVAL  OF  ITEMS  FROM  LIST  IS  ACCOMPLISHED  AS  PART  OF  STEPS  3 
AND  4  AT  STATEMENT  248  BELOW.) 

200  IF(NITEMS.EO.O)  RETURN 
LR  *  LOCITMlli 

READ  ( 1 ' LR  )  (INFO(L)»L*lfLINFO) 

NITEMS  =  N1 TEMS  -  1 

IF  INI TEMS.EQ.O)  GO  TO  225 
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C  PUSH  DOWN  LIST. 

00  210  L  -  l.NITEMS 
ZLINIT(L)  >  ZLINITCL+l) 

210  LOCITM(L)  «  LOC I TM( L+l ) 

LOCITMINITEMSM  )  «  LR 

C  SAVE  OATA  FOR  STEP  4. 

22S  XBKS  *  XBK 

OS  *  OIJBK) 

ZBKS  *  ZBK 
CS  *  C(JBK) 

ZPS  *  OIJBK) 

OSUHS  «  OSUM 
JBKS  «  JBK 

C  STEP  3  INITIALIZATION. 

NCASE  «  1 
QCJBKI  «  XBK 

C(JBK)  *  IZBK-ZPS)/(XBK-P( JBK) ) 

GO  TO  240 

C  STEP  4  INITIALIZATION. 

235  NCASE  «  2 

P(JBKS)  *  XBKS 
O(JBKS)  ■  OS 
Z 0  *  PHI ( JBKS (OS ) 

C ( JBKS )  «  ( ZQ-ZBKS ) / ( QS-XBKS ) 

DC JBKS)  *  ZBKS 
OSUM  *  OSUM  -  ZPS  ♦  ZBKS 
240  NITER  «  NITER  ♦  1 

C  FIND  OPTIMUM  SOLUTION  ASSUMING  LINEARIZED  COSTS. 

300  00  305  J  -  1,N 
305  XLINtJt  -  0. 

I  «  0 

L  *  l 

310  I  ■  I  ♦  l 

ZI  «  FINF 
H  ■  0 

L  »  L  ♦  1 

320  H  »  H  ♦  i 

ZIH  -  0. 

LIHl  *  L  ♦  l 
330  L  -  L  ♦  1 
JJ  *  JIN(L) 

IF ( JJ.LE.O)  GO  TO  340 
ZIH  «  ZIH  ♦  X  I N ( L ) *C I J J ) 

GO  TO  330 

340  IFIZIH.GE.Zl )  GO  TO  345 
Zl  «  ZIH 
HI  ■  H 
L 1 1  «  LIHl 
LI2  *  L  -  1 
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345  lMJJ.60.0i  GO  TO  320 
HILIN(I)  «  H! 

00  350  L350  »  LIltLI2 
J350  *  J 1 N( L350 ) 

350  XL  I N ( J350 )  *  XLIN( J350)  ♦  XfN(L350) 

IF(JJ.EO.-l)  GO  TO  310 

C  COMPUTE  LINEAR  COST  ZLIN  AND  TRUE  COST  ZMIX  OF  THIS  SOLUTION. 

ZLIN  *  DSUM 
ZMIX  *  0. 

00  241  J  *  1 « N 

ZLIN  «  ZLIN  ♦  (XLIN(J)-P(J)  ) *C( J ) 

Z( J)  *  PHI (JyXLINlJ) ) 

241  ZMIX  =  ZMIX  ♦  Z(J) 

IF(ZMIX.GE.ZBEST)  GO  TO  250 

C  IF  THIS  IS  A  BETTER  SOLUTION  THAT  THE  BEST  SO  FAR,  REPLACE  IT. 

ZBEST  =  ZMIX 
00  245  I  «  1 , M 
245  HIBFSTI I J  *  HILINI I ) 

DO  247  J  *  l,N 

247  XBEST(J)  =  XLIN(J) 

C  CHECK  LIST  ZLINII  FOR  REMOVABLE  ITEMS. 

248  IF(NITEMS.EO.O)  GO  TO  250 
IFIZLINITINITEMSI.LT. OMEPS*Z BE ST)  GO  TO  250 
NITEMS  *  NITEMS  -  l 

GO  TO  248 

C  DETERMINE  IF  THIS  CASE  SHOULD  BE  PUT  ON  THE  LIST  OR  IGNORED. 

250  IFIZLIN.GE .ZBEST 40MEPS )  GO  TO  260 
C  DETERMINE  INTERVAL  TO  BE  8R0KEN  FOR  THIS  CASE. 

EMAX  *  0. 

DO  255  J  *  l,N 

E  *  Z(J)  -  Of  J )  -  C(  J)MXLIN(J)-PCJ)  ) 

IF (6. LE. EMAX)  GO  TO  255 

EMAX  -  E 

JBK  «  J 

XBK  *  XLIN(J) 

ZBK  *  ZIJ) 

255  CONTINUE 

I  F ( EMAX.EQ.O.  )  GO  TO  1100 

IF ( XBK.LE.P! JBK) ♦( 1 . +XE PS ) .OR. XBK .GE . Q ( JBK ) * (  l.-XEPS) )  GO  TO  1100 
C  SORT  THE  VALUE  ZLIN  INTO  THE  LIST  ZL1NIT,  SAVE  THIS  CASE  ON  OISK. 

IFINITEMS.GE.ITEMAX)  GO  TO  1000 
NITEMS  «  NITEMS  ♦  1 

IFINITEMS.GT.MAXITM)  MAX  I TM  =  NITEMS 
LW  *  LOCITMINI TEMS) 

IT  *  NITtMS 

256  IF! IT.LE.l)  GO  TO  257 

I F  C  ZL  I N I  T ( IT-1). LE. ZLIN)  GO  TO  257 
ZLINITUI)  •  ZLINITI  IT-1  ) 

LOCITMI IT)  =  LOCITMI IT-1  ) 
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IT  *  IT  -  1 

GO  TO  256 

257  ZllNITUn  ■  ZL  I  N 
LOCI TM( I T )  «  LW 

WRITE  (  I  •  L  W I  (  INFO(L),L*l»LINFO) 

260  IFINCASE.EQ.2)  GO  TO  200 
GO  TO  235 

900  WRITE  (6,901)  M,N,L 

901  FORMAT ( '  SUBROUTINE  CAVAP  ENTERED  WITH  IMPROPER  DATA. '/lOX'M  «'I4 

1  *,  N  *'14',LAST  LOCATION  IN  INPUT  LISTS  JIN  AND  XIN  CHECKED  «' 

2  15) 

RETURN  1 

1000  WRITE  (6,1001)  N1TEMS, NITER 

1001  FORMAT ( *  SUBROUTINE  CAVAP  DISK  STORAGE  OF  *15*  ITEMS  EXHAUSTED  AFT 
lfcR'I6'  ITERATIONS.') 

RETURN  2 

1100  WRITE  (6,1101)  JOK, NITER 

1101  FORMAT! •  SUBROUTINE  CAVAP.  ERROR  DETECTED  IN  COMPONE NT ' 1 3 •  AFTER' 
115'  ITERATIONS.'/'  MAY  BE  CAUSED  BY  ROUND-OFF  ERRORS  OR  NONCON 
2CAVITY  OF  FUNCTION  PHI.*) 

RETURN  2 
END 

TOTAL  MEMORY  REQUIREMENTS  002704  BYTES 
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